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In neutral graphene dots the Fermi level coincides with the Dirac points. We have investigated 
in the presence of a magnetic field several unusual properties of single electron states near the 
Fermi level of such a rectangular-shaped graphene dot with two zigzag and two armchair edges. 
We find that a quasi- degenerate level forms near zero energy and the number of states in this level 
can be tuned by the magnetic field. The wavefunctions of states in this level are all peaked on 
the zigzag edges with or without some weight inside the dot. Some of these states are magnetic 
field-independent surface states while the others are field-dependent. We have found a scaling result 
from which the number of magnetic field-dependent states of large dots can be inferred from those 
of smaller dots. 

PACS numbers: 



INTRODUCTION 




FIG. 1: Finite graphene layer with zigzag and armchair edges. 
There are equal number of A and B carbon atoms. The 
graphene layer has reflection symmetries about horizontal and 
vertical lines that go through the center of the layer. A mag- 
netic field is present perpendicular to the layer. 



Graphene dots have a great potential for many ap- 
plications since they are the elemental blocks to con- 
struct graphene-based nano devices. It is possible to cut 
graphene sheet [1] in the desired shape and size [2], and use 
it to make quantum dot devices. In such devices it may 
be possible to realize experimentally zigzag or armchair 
boundaries. 

Graphene systems possess several unusual physical 
properties associated with the presence of the Dirac 
points. For example, compared to ordinary Landau lev- 
els of quasi-two-dimensional semiconductors the lowest 
Landau level (LLL) of graphene is peculiar since it has 
zero energy that is independent of magnetic field (sj-Q- 
Moreover, wavefunctions of the LLL are chiral, i.e., the 
probability amplitude of find the electron on one type 
carbon atoms is zero. There are other graphene systems 
with zero energy states. Semi-infinite [8[ or nanoribbon 



graphene [iiliol with zigzag edges along the x-axis develop 
a flat band of zero energy chiral states. These states are 
surface states and are localized states at the zigzag edges 
with various localization lengths [8|. The zigzag edge and 
the LLL states have zero energy because their wavefunc- 
tions are chiral. Effects of a magnetic field on graphene 
Hall bars have been investigated recently, and some zero 
energy chiral states are found to be strongly localized on 
the zigzag edges in addition to the usual LLL states [n|- 

13. 

One may expect that the degeneracy of chiral states 
with zero energy will be split when quantum confine- 
ment effect is introduced in a graphene dot p!5l-[l7| . How- 
ever, the splitting of these energies may be unusual 
in some graphene dots^-lHI- Recently the magnetic 
field dependence of these levels in a gated graphene 
dot was investigated experimentally [21] . Effects of var- 
ious types of edges have been also investigated: zigzag- 
edged dots, armchair-edged dotsfiol, |22|-[2J, and rectan- 
gular graphene dots with two zigzag and two armchair 
edges j25l [26] have been studied. Armchair edges cou- 
ple states near K and points of the first Brillouin 
zone and generate several mixed chiral zigzag edge states 
with nearly zero energies. In the rectangular dots the 
number of these states. A/"/, may be determined from the 
condition that the x-component of wavevectors satisfies 



1/Ly < kx^n < 7r/3a, where kx,n = 17 " fj^ ^ = 



ao 



is the length of the unit cell, and n = 0, ±1, ±2, ... (the 
nearest neighbor carbon-carbon distance is ao = 1.42 A, 
the horizontal length of the dot is Lx = V^Mao with 
M number of hexagons along the x-axis, and the vertical 
length is Ly = ao{3N-\-2) with N the number of hexagons 
and {N + 1) carbon bonds along the y-axis. See Fig(T]). 
This condition implies that the integer n is given by 



{3N + 2)7r 



2M 



<n< M. 



(1) 
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The effective mass approximation wavefunctions of these 
surface states are derived in Ref.[26]. 

We investigate how properties of rectangular dots 
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FIG. 2: Number of nearly zero energy states that are induced 
by a magnetic field follow a scaling curve when plotted as a 
function of (^^^jJ^^2/x/3)0 different rectangular shaped 
graphene dots. Here 5 = O.OleV. 



change in the presence of a magnetic field in the regime 
where Hofstadter-butterfly effects |l8|, UH, [111 are negligi- 
ble. In neutral graphene dots the Fermi level has zero 
energy, and, consequently, magnetic, optical, and STM 
properties are expected to depend on the number of avail- 
able states near zero energy. Our investigation shows 
that the wavefunctions of nearly zero energy states are 
all peaked on the zigzag edges with or without appre- 
ciable weight inside the dot. This may be understood 
as mixing of LLL and surface states by armchair edges 
of the square dot through intervalley scattering, which 
is unique to the square dot (This will be explained in 
Sec.IV). 

Our study shows that the number of states within the 
energy interval S around zero energy is given by 



Nt{^)=Ni^NdW 



(2) 



(the energy S is typically less than the quantization en- 
ergy of a rectangular dot, which can be estimated using 
the Dirac equation: jkmin ^ ^r^? where 7 = >/3ta/2 
with the hopping energy t). Here magnetic edge states 
are not included since their energies are larger than S. 
N]j{(j)) is the number states at zero magnetic field that 
merge into the energy interval 6 around zero energy as 
the dimensionless magnetic flux cj) = increases. This 
effect provides a means to control the number of states at 
the Fermi level. There are other states with nearly zero 
energies at (/) = 0, which remain so even at 7^ 0. The 
localization lengths of these states are shorter than the 
magnetic length. We denote the number of these states 
hy Ni. 

Our numerical results indicate that for relatively small 
rectangular-shaped graphene dots with size less than of 
order 10^ A the number of states at the Fermi level display 
a negligible magnetic field dependence for values that are 
usually accessible experimentally {B < lOT corresponds, 
in dimensionless magnetic flux, to = < 10~^, where 
<l>o = hc/e and ^ = BAh with the area of a hexagon 



Ah 



3^/3^2 



0; 



On the other hand, in larger dots this 



dependence is significant. However, it is computationally 
difficult to investigate large rectangular-shaped graphene 



dots since the number of carbon atoms increases rapidly 
with the size. We have found a scaling result from which 
one can infer results for larger dots from those of smaller 
dots. For different rectangular-shaped graphene dots and 
values of (j) it can be described well by the following di- 
mensionless form 



(3) 



where i 



ao{L, 

is the magnetic length and f{x) is a 



(47r0)i/2 

scaling function, see FiglS The total number of hexagons 
in the dot is = L^Ly/Ah. Our numerical result shows 
that the dependence of N]j{(j)) on (j) is initially non-linear 
in the regime where the diameter of the cyclotron motion 
is comparable to the system length, 2£ ^ Lx^y. 



II. NUMBER OF STATES IN THE 
QUASI-DEGENERATE LEVEL 



Our Hamiltonian is 

H = — tijc\cj^ 

<i,j> 



(4) 



where t 



if jy A-dr 

te ^0 are the hopping parameters and 



c- creates an electron at site i. 



Here we use a Lan- 
dau gauge A = 5(— y,0,0). The summation < z,j > 
is over nearest neighbor sites and t = 2.5eV. The eigen- 
state with eigenenergy Cn is denoted by (j)e^{R)^ where R 
labels each lattice point. Because of electron- hole sym- 
metry eigenvalues appear in pairs of positive and nega- 
tive values (e, — e), and the probability wavefunctions of a 

pair of states, (|^e^(^)P, \^-eri{^)\^)i identical. Our 
numerical results are consistent with this. 

FigslSja) and (b) display the energy spectra near zero 
energy for (/) = and 0.01. At = there are approx- 
imately 20 states within |en| < O.OleV, consistent with 
the analytical result of Eq.([T]). At ^ = 0.01 the numerical 
value is increased to 24. FigHJ^a) shows how some energy 
levels en at (/) = change as a function of the magnetic 
flux (j). These energy levels do not anticross. We ob- 
serve that nearly zero energies at (/> = do not change 
noticeably in magnitude as (j) varies. There are Ni such 
localized surface states. On the other hand, we find that 
as (j) increases non-zero energies become smaller and move 
closer to zero. This implies that, for a given energy inter- 
val the number of states in it, A^t(^), increases with (j). 
From FigHJb) we see that it displays a non-linear depen- 
dence on (j). For a large dot of size 50 x 50nm^ a similar 
dependence of Nt on B is seen, as shown in Fig|5l Non- 
linear dependence occurs in the regime 2^/La^ ~ 0.5. As 
a test of our numerical procedures we have verified that 
the sum of A^d(0) = Nt{(I)) — Ni and the number of mag- 
netic edge states is equal to the total bulk Landau level 
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FIG. 5: Results for a dot with size 50 x 50nm^ . Dependence 
of Nt/2 on B for ^ = O.OleV. 



FIG. 3: (a) Eigenenergies Cn at = 0. Quasi-degenerate 
states are present near zero energy. Size of the dot is 74 x 
71 A^. A quantization energy of order jkmin ^ 0.03eV can be 
seen as an excitation gap near zero energy, (b) Eigenenergies 
en at 6 = 0.01. 
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FIG. 4: (a) Some energy levels Cn change as increases while 
some do not. (b) The total number of states within the energy 
interval S around zero energy at a finite value of (j). Size of 
the dot is 74 x 71 



The ratio between the number of nearly zero energy 
states induced by the magnetic field and the number 
of hexagons, should depend on a dimensionless 

quantity consisting of a combination of ^, ao, and 
Ly, which are the important parameters of rectangu- 
lar graphene dots. The lengths Lx and Ly should ap- 
pear as Lx + Ly so that for rectangular-shaped graphene 
sheets No/Nh remains the same when Lx and Ly are ex- 
changed. These considerations lead us to the dimension- 
less variable — r/ \ = t-tttj — t^t^^tw^, see Eq.JSl). 
We are especially interested in the regime where the di- 
ameter of the cyclotron orbit is comparable to the sys- 




::S::::::::::::::::::::::::!:::::::ct:::o:::::::::::::: 



. v.v.v^.-;:" • v.;-; 
.^. /Xv/X* •■•Xv.-X • 



FIG. 6: Size of the dot is 74 x 71 (a) The probability 
wavefunction of the state with eio27 = — 0.07eV at 0=0. The 
length unit is a. (b) Profile of z-component of pseudospin: 
sizes of red (blue) dots represent probabilities of occupying A 
(B) carbon atoms. Note that blue (red) dots are dominant 
in the upper (lower) of dot. When the probabilities are less 
than 0.00001, the radius of the dots is set to the smallest 
value. The upper and lower horizontal edges represent zigzag 
edges, (c) The probability wavefunction of the state with 
61027 = -2.0 X 10-^eV at i/L^^ = 0.12 (0 = 0.01). (d) Profile 
of z-component of pseudospin for n = 1027 at = 0.01. 



tem length 2£ ^ Lx^y Note that in this regime many 
cyclotron orbits get affected by the presence of the edges 
and corners of the rectangular dot. Since we must also 
assume that Hofstadter effect is negligible the validity 
regime of Eq.([3|) is ao <C 2^ < Lx^y Note also that the 
scaling function f{x) should be different for each 5. 



III. WAVEFUNCTIONS OF 
QUASI-DEGENERATE STATES IN A 
MAGNETIC FIELD 

We first show how the wavefunction of a non-zero en- 
ergy state at = changes into a state with nearly zero 
energy as (j) increases. Consider the probability wave- 
function for n = 1027 at a finite (j) = 0.01, as shown in 
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FIG. 7: Size of the dot is 74 x 71 (a) The probabihty 
wavefunction of the state with eio45 = 5.38 X 10"^eV at = 
0. (b) Profile of z-component of pseudospin for n = 1045 
at = 0. (c) The probability wavefunction for n = 1045 
at i/La^ = 0.12 (0 = 0.01). (d) Profile of z-component of 
pseudospin for n = 1045 at = 0.01. 




FIG. 9: (a) The probability wavefunction of the state with 
£1053 = 0.144eV at i/L^, = 0.12 (0 = 0.01). Size of the dot 
is 74 X 71 A^. (b) Profile of z-component of pseudospin of the 
same state. 
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FIG. 8: Size of the dot is 7Ax7lX^ and l/L^ = 0.09 (0 = 0.02) 
(a) The probability wavefunction of the state with eio27 = 
—2.04 X 10~^eV. (b) Profile of z-component of pseudospin of 
the same state. These results should compared with those in 

Figil 



FigEl^c). It is localized on the zigzag edges with a fi- 
nite probability inside the dot. On the armchair edges 
the wavefunction is vanishingly small. The wavefunc- 
tion has changed significantly from the = result, 
see Figin^a), and also its energy has changed from - 
O.OTeV to —2.0 x 10~^eV. The values of the z-component 
of the pseudospin, FiglGfb) and FiglGjd), are larger on 
the zigzag edges at (j) = 0.01 compared to the result at 
= 0. The probability wavefunction of another state 
with nearly zero energy is shown in FigH^c) at = 0.01. 
We see that the result is somewhat different from the zero 
field result of Figd^a), which displays a localized state 
with the localization length comparable to the unit cell 
length a. Now there is a finite probability to find an elec- 
tron inside the dot while the probabilities on the zigzag 
edges are reduced. Note that the energy of this state has 



changed from 5.38 x 10~^eV to 4.5 x 10~^eV when is 
changed to 0.01 from zero. The pseudospin profiles are 
shown in FigsITfb) and (d). Fig 15] shows a probability 
wavefunction at a smaller value of (^/L^^y = 0.09 (corre- 
sponding to = 0.02), and we see that the wavefunction 
is less localized on the zigzag edges and LLL character is 
more pronounced in comparison to result of ^/ = 0.12 
(FigEl^c)). All N]j{<p) states have similar properties men- 
tioned above with finite probabilities of finding an elec- 
tron inside the dot. There are also N£ zigzag edges states 
that are more strongly localized on the edges with local- 
ization lengths comparable to a. When an electron is in 
one of these states the probability of find the electron 
away from the edges is practically zero. The energies of 
these states are less then 10~^^eV. We can summarize 
our results as follows: all nearly zero energy states are 
localized on the zigzag edges with or without some weight 
inside the dot. 

When the magnetic length is much smaller than the 
system size magnetic edge states can be formed, see FiglHl 
The probability wavefunction of a magnetic edge state 
with 61053 = 0.144eV is shown in Fig|9](a). It is a mixture 
of ordinary magnetic and zigzag edge states. The proba- 
bility wavefunction decays from the armchair edges while 
it is strongly peaked on the zigzag edges. Pseudospin 
expectation values on each zigzag edge display opposite 
chiral behavior. On armchair edges the chiralities are 
more or less evenly mixed, see Fig|9jb). 



IV. DISCUSSIONS AND CONCLUSIONS 

We now explain qualitatively how mixed states of 
Figin^c) and Fig|7l(c) can arise. An infinitely long 
zigzag nanoribbon in a magnetic field has nearly zero 
energy surface states that are localized on the edges 
in addition to ordinary lowest Landau level states, see 
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FIG. 10: Cross section of probability wavefunctions of a 
nanoribbon with infinitely long zigzag edges along the x-axis 
in the presence of a perpendicular magnetic field (a Landau 
gauge is used). represents a localized surface state. Two 

examples of LLL states |'02|^ and I'^sl^ are also shown. These 
states all have nearly zero energies. 



FigJTOl The properties of these states are given in 
Refs.ji, [O, EM HI: LLL states of valley K (K^ are 
of B (A) type and localized surface states have a mixed 
character between A and B. The surface states can have 
various localization lengths but the minimum value is 
of order the carbon-carbon distance an 1121. The arm- 
chair edges couple K and vallevsfpl, I 111 |26|, and, con- 
sequently, surface and LLL states of a nanoribbon can 
be coupled and give rise to mixed states with significant 
weight on the zigzag edges and inside the dot, as shown 
in FigEfc) and Fig|7fc). In addition, these mixed states 
should display a significant occupation of both A and B 



carbon atoms inside the dot since LLL states of different 
chiralities are coupled by the armchair edges. Our nu- 
merical result is indeed consistent with this expectation, 
see Figini (d). As the ratio i/Lx,y takes smaller values 
the nature of these states become more like that of LLL 
states (see FigEj). 

We have investigated, in the presence of a mag- 
netic field, quasi-degenerate states of rectangular-shaped 
graphene dots near the Dirac points. Some of these states 
are magnetic field independent surface states while the 
other states are field dependent. We find numerically 
that the wavefunctions of these states are all peaked near 
the zigzag edges with or without significant weight inside 
the dot. The physical origin of the presence of a signifi- 
cant weight is the coupling between K and K' valleys due 
to the armchair edges. This effect is expected to survive 
small deviations from perfect armchair edges as long as 
they provide coupling between different valleys. Experi- 
mentally the dependence of N]j on (j) may be studied by 
measuring STM properties [l^ or the optical absorption 
spectrum as a function of magnetic field [28- 30]. In fab- 
ricating rectangular dots a special attention should be 
given to the direction of armchair edges since the prop- 
erties of dot may depend on itfss']. 

Acknowledgments 

This work was supported by the Korea Research Foun- 
dation Grant funded by the Korean Government (KRF- 
2009-0074470). In addition this work was supported by 
the Second Brain Korea 21 Project. 



[1] For recent reviews see: T. Ando, J. Phys. Soc. Jpn. 74, 

777 (2005); A. K. Geim and A. H. MacDonald, Phys. 

Today 60 35 (2007); A. H. Castro Neto, F. Guinea, N. 

M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. 

Phys. 81, 109 (2009). 
[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004); H. Hiura, Appl. Surf. 

Sci. 222, 374(2004); A. K. Geim and K. S. Novoselov, 

Nat. Mater. 6, 183 (2007). 
[3] Y. Zhang, Y.W. Tan, H.L. Stormer, and P. Kim, Nature 

(London) 438, 201 (2005). 
[4] K. S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, 

M.I. Katsnelson, I.V. Grigorieva, S.V. Dubonos, and 

A. A. Firsov, Nature (London) 438, 197 (2005). 
[5] J.W. McClure, Phys. Rev. 104, 666 (1956). 
[6] Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
[7] S. G. Sharapov, V. P. Gusynin, and H. Beck, Phys. Rev. 

B 69, 075104 (2004). 
[8] M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusak- 

abe, J. Phys. Soc. Jpn. 65 ,1920 (1996). 
[9] L. Brey and H. A. Fertig, Phys. Rev. B 73, 195408 (2006). 
[10] Y. W. Son, M. L. Cohen, and S. G. Louie, Phys. Re v. Lett. 

97, 216803 (2006). 
[11] L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 
[12] A. H. Castro Neto, F. Guinea, and N. M. R. Peres, Phys. 



Rev. B 73, 205408 (2006). 
[13] V. P. Gusynin, V. A. Miransky, S. G. Sharapov, 1. A. 
Shovkovy and C. M. Wyenberg, Phys. Rev. B 79, 115431 
(2009). 

[14] M. Arikawa, Y. Hatsugai, and H. Aoki, Phys. Rev. B 78, 
205401 (2008). 

[15] A. DeMartino, L. DellAnna, and R. Egger, Phys. Rev. 

Lett. 98, 066802 (2007). 
[16] G. Giavaras, P. A. Maksym, and M. Roy, J. Phys.: Con- 
dens. Matter 21, 102201 (2009). 
[17] W. Hausler, and R. Egger, arXiv:0905.3667^1. 
[18] S. Schnez, K. Ensslin, M. Sigrist, and T. Ihn, Phys. Rev. 

B 78, 195427 (2008). 
[19] Z. Z. Zhang, K. Chang, and F. M. Peeters, Phys. Rev. B 

77, 235411 (2008). 
[20] P. Recher, J. Nilsson, G. Burkard, and B. Trauzettel, 

Phys. Rev. B 79, 085407 (2009). 
[21] J. Giittinger, C. Stampfer, F. Libisch, T. Frey, J. 

Burgdorfer, T.Ihn, and K. Ensslin, cond-mat arXiv: 

0904.3506v2. 

[22] J. Fernandez-Rossier and J. J. Palacios, Phys. Rev. Lett. 

99, 177204 (2007). 
[23] M. Ezawa, Phys. Rev. B 76, 245415 (2007). 
[24] O. Hod, V. Barone, and G. E. Scuseria, Phys. Rev. B 77, 

035411 (2008). 

[25] B. Trauzettel, D. V. Bulaev, D. Loss and G. Burkard, 



6 



Nature Phys. 3, 192 (2007). 
[26] C. Tang, W. Yan, Y. Zheng, G. Li, and L. Li, Nanotech- 

nology, 19, 435401 (2008). 
[27] Y. Niimi, H. Kambara, T. Matsui, D. Yoshioka, and H. 

Fukuyama, Phys. Rev. Lett. 97, 236804 (2006). 
[28] M. O. Goerbig, J. -N. Fuchs, K. Kechedzhi, and V. L 

Fal'ko, Phys. Rev. Lett. 99, 087402 (2007). 
[29] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. 



Rev. Lett. 98, 157402 (2007). 
[30] H. Hsu and L. E. Reichl, Phys. Rev. B 76, 045418 (2007). 
[31] M. Kohmoto and Y. Hasegawa, Phys. Rev. B 76, 205402 

(2007). 

[32] L.Brey and H.A. Fertig, Phys. Rev. B 75, 125434 (2007). 
[33] A. R. Akhmerov and C.W.J. Beenakker, Phys. Rev. B 
77, 085423 (2008). 



